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INTRODUCTION 


Volume six in the college physics series has these prograas: 


THERMODYNAMIC PROCESSES CALORIMETRY 
THERMODYNAMIC CYCLES IDEAL GAS CYCLES 
HEAT ENGINES - THEORY MOLECULAR MOTION AND PRESSURE 


HEAT ENGINES - APPLICATIONS DEMUFF IN 
plus several other text and binary files used by “Ideal Gas Cycles". 


IDEAL GAS CYCLES was written by Paul Stephenson. MOLECULAR MOTION 
AND PRESSURE was written by Mark Cross. Steve Kama wrote all of the 
others, 


This diskette uses Apple DOS 3.3. If you have an older disk drive 
with DOS 3.2 then you will have to BRUN DEMUFFIN to convert the whole 
disk. Jake a blank diskette initiallized under DOS 3.2, plus this 
physics programs diskette, to a friend who has DOS 3.3. He will show 
you how to convert the prograas to D0S 3.2. 
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Isobaric process on a V-T diagram 
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THERMODYNAMIC PROCESSES 


The program draws isobaric, adiabatic, isometric, and isothergal 
Processes and asks questions about them. It begins with a true-false 
quiz about the definitions of these processes. Next it draws 
families of curves for an ideal gas with different nushers of 
molecules going through each process. The curves are repeated on 
P-¥, P-T, and V-T graphs. Examples are shown below. Finally it asks 
the viewer three questions like this one: | 


UNDER A GRAPH OF A SLANTED LINE ON A V-T DIAGRAM ) 
IS THIS: 1. ISOTHERMAL ON V-T DIAGRAM 

2. TSDBARIC ON V-T DIAGRAM 

3, ISORARIC ON P-T DIABRAM 

ANSWER 1, 2, OR 3: 


Isothermal process on a P-V graph 
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THERMODYNAMIC CYCLES 


This program assumes the viewer has seen "Theraodynasic 
Processes". It builds on that program to combine the basic processes 
into cycles. 


The first cycle illustrated has isobaric expansion, then 


isometric cooling, then isoharic compression, and finally isometric 
heating. It is a rectangle on a P-¥ diagram. The program shades in 
the area below the isobaric expansion and then it subtracts the area 
under the compression to get the net work. 
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The program studies a second cycle gade by an isobaric expansion 
combined with an isometric cooling followed by isothermal 
compression. It draws the cycle as before. Then it discuss finding 
the net work by integrating PdV with the ideal gas equation. 


The program briefly mentions the Carnot cycle. It says that the 
cycle is made of 2 isothermal processes and 2 adiabatic processes, 
and that it results in the maxisum possible efficiency. 


HEAT ENGINES — THEORY 


This program is meant to be used for further study after the 
Carnot cycle has been presented in a physics lecture. It begins with 
a reminder that the Carnot cycle gives the best possible efficiency 
of a heat engine or puap, but real engines have less than this 
efficiency. 


The program shows an animated heat flow diagram for an engine 
with the hot reservoir on the top, the cold on the bottom, and work 
being extracted at the side. Then it says that a heat puap is the 
opposite, and draws a pump with heat moving from the cold to the hot 
area, 


Using the standard notation with “1" being the hot temperature, 
heat engines have efficiency = W/@1. Heat pumps have coefficient af 


performance kK = Q2/W. These formulas are also presented in sentence 
form, Next come three T/F questions: 


1, Use efficiency for heat pumps? T or F? 
2, @2 is apasured in degrees Kelvin? T or F? 
3, kK =W / @? T or F? 
The viewer i5 expected to get all three answers right. 
These final two formulas are presented with some text: 
W= O1 - O2: and @1/82 = T1/T2. Then the progras uses them to solve 
this problem: 
A heat pump operates between -!0 deorees and 25 degrees 


celsius. How much work must be done to extract 200 calories of 
thermal energy from the low teaperature sink? 


HEAT ENGINES — APPLICATIONS 


This program applies the theory of heat engines to particular 
probleas involving ideal engines and refrigerators. The viewer 


should have a pencil, paper, and calculator. It starts with this 


review and summary: 


HIGH TEMP SOURCE (KELVIN) 
LOW TEMP SINK (KELVIN) 

HEAT FLOW FROM HIGH TEMP 

Q2 = HEAT FLOW FROM LOW TEMP 

W = WORK DONE BY THE ENGINE 

EFF= EFFICIENCY = W / Qi 

K = HEAT PUMP PERFORMANCE = Q2 / W 


RELATIONSHIPS: W=01- 02: @1 / @2= Ti / 72 


PROBLEMS SOLVED: 


(1) How auch useful work can be extracted from an ideal heat engine 
operating between 200 degrees and 105 degrees celsius if 1 kcal 
of energy is supplied? 


(2) An ideal refrigerator extracts heat from a storage chamber at 
-15 degrees and exhausts it at 40 degrees celsius. How many 
kcal are extracted if 100 joules of work are supplied? 


(3) What is the aaxiaue possible efficiency of a steam engine if the 
teaperature of the input steaw is 175 degrees and the exhaust is 


7a celsius? 


(4) fn engine absorbs 500 kcal of heat at 600 degrees kelvin and 
exhausts heat at 300 kelvin. What is the efficiency and how 
much heat is exhausted, assuming an ideal engine? 


(5) A Carnot engine operates between 20 and 500 degrees celsius, 
How much energy aust be supplied to produce 10,000 joules af 


work? 
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CALORIMETRY 


When you start with a hot sample and put it into a cup of cool 
liquid, both cose to the same teaperature. The progras shows this process 
with pictures, 


The law of conservation of energy says that the heat lost by the hot 
object equals the heat gained by the cold cup and the liquid. Heat energy 
for these problems is calculated froa specific heat times aass times 
temperature change. After presenting this theory, the program asks the 
viewer to make a blank chart like the one below. Some teachers aight want 
to leave blank copies of this chart beside the coaputer, 
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Next the progras solves the six probleas listed below. After the 
third problem, the computer teaches the concept of the water equivalent of 
an object. It uses this idea in the fourth problem. It teaches about the 
heat of fusion between probleas 4 and 5, 


Problea 1, A 600 ga copper cup (c = .093) holds 1500 ga of water at 20 
degrees C. A 100 ge sample of iron (c = .11) at 120 degrees Cis put in 
the water. What is the final teaperature? 


Problem 2, A 300 ga glass container (c = .14) holds 500 ga of water at 
30 degrees C. When 500 ga of alcohol at 75 degrees are added, the sixture 
reaches a final teaperature of 44 degrees C. What is the specific heat of 
the alcohol? 


Problea 3, A 450 ga piece of lead is heated to 100 degrees C and then 
dropped into a 500 gm copper {c = .093) calorimeter cup which contains 100 
qa of water at 10 degrees C. After @ixing, the teaperature of the water 
has increased 11.1 degrees. Find the specific heat of the lead. Hint: the 
equilibrium temperature is not {1.1 degrees! 


Problea 4. A calorimeter (water equivalent = .03 kg) holds 500 ga of 
water at 15 degrees C. A 540 gm saaple of copper {c = .093) is added. If 
the final temperature of the water is 22,5 degrees (, find the initial 
temperature of the copper, 


Problema 5. A 30 gm piece of ice at zero degrees Cis placed in 210 ga 
of water at 30 degrees C. Neglecting heat losses to the environaent, what 
is the final teaperature? 


Problem 4. 500 ga of lead {c = .03) pellets at 100 degrees C are poured 


in a hole in a large piece of ice. How auch ice is aelted before the lead 
cools to zero? 
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IDEAL GAS CYCLES 


This progras simulates an ideal gas enclosed in a cylinder 
with a moveable piston. Several processes may be perforaed on a 
monatomic, diatomic, or polyatomic ideal gas. Although the use of 
a textbook by the student is assumed, he can sSumson some 
definitions and general instructions by pressing the "ESC" key when 
the computer is waiting for a nugerical response, but not from the 
@ain menu, 


The aenu gives seven choices: 
{. ISOTHERMAL EXPANSION 3. RECTANGULAR CYCLE 
2. ADIABATIC EXPANSION & CARNOT CYCLE 
3. ISOBARIC EXPANSION 7, REFRIGERATOR CYCLE 


4, TSOCHORIC PROCESS 8. QUIT 


The user picks 


Choices 1-4 allow the user to pick the beginning and end 
points of the process. Pressure is in atgospheres and volume in 
liters, making the Ideal Gas equation: P Y = .0823 7. The program 
plots the chosen process on a graph like the one below. The gas 
expands and contracts in the cylinder at the right side of the 
picture and the heat reservoir is connected or disconnected as 
required. The user can plot several processes on one graph to 
build up a whole cycle. 
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Choices 3-/ allow the user to specify the four quantities that 
determine a closed cycle. For the Carnot and refrigerator cycles, 
the user picks the upper temperature and three of the volumes. 
Then the computer goes through the cycle, soving the gas up and 
down in the cylinder and calculating the heat and work as it goes. 
At the end it prints the efficiency of the whole cycle. 


Ti, VA, VB, Vc 


MOLECULAR MOTTON AND GAS PRESSURE 


This program simulates the motion of ideal gas solecules. 
They bounce around at random in a cylinder, sometiaes hitting the 
Pressure gauge on the right side. Every hit on the qauge causes a 
temporary small pressure. The pressure fluctuations are saoothed 
out a little and then graphed in the middle of the screen. 


This deaonstration Prepares students for the derivation of the 
relation between gas kinetic theory and the average pressure. They 
can watch wall collisions and their effect in making pressure 
fluctuations, 


The program is easier to use if people know about the repeat 
key. It can be held down together with another key such as "M" for 
acre molecules. 


The gas particles have a Maxwellian distribution of speeds, 
The whole theoretical distribution is divided into twenty bins of 
equal probability. That is, the area in each bin is the same on 
the graph of n{v) versus v. New molecules are assigned to a randos 
bin to get their speed. They are then placed in the center of the 
cylinder and started moving in the general direction of the lower 
right corner. 


NI 


The temperature originally is "proportional ta the average 
speed squared. To save time, later temperatures are found froe 
changes in the rene If the speed of every aolecule increases by 
a factor of 1.2 then the teaperature goes up by the square of this, 
or 1.44 times. At the same tiae the whole Maxwellian distribution 
in the graph above would be aultiplied by 1.2. 
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PROGRAMMING DETAILS 


THERMODYNAMIC PROCESSES and THERMODYNAMIC CYCLES use text and 
the high resolution screen HGR. They have no unusual programming 
tricks. Their aain program flows from line 100 through to about 
800 with very few branches out to subroutines, except for frequent 
juaps to line 10050 to wait for the user’s carriage return, 


ENGINES: THEORY uses text and animation on the low resolution 
screen. ENGINES: APPLICATIONS is straight text problea oon with 
na graphics animation, 


CALORIMETRY uses some low resolution graphics plus a lat of 
text. Its structure is the same as the above four prograas; a ain 
program that flows straight through from line 100 to about 800, 


IDEAL GAS CYCLES is extremely coaplex. It uses both high 
resolution screens. The Basic code resides above the second screen 
HGR2 when it runs. Machine language routines and other smal! itess 
are stored in the normal program space between aesory locations 
2048 and 8192. Al] this configuration is set up by the disk file 
EXEC AUTO. This file loads some binary files, sets up the start of 
program pointer, and then RUNs the main program: GAS.COM.4. Once 
IDEAL GAS CYCLES is running, it gives beautiful high resolution 
graphics using the character generator in Apple’s DOS TOOL KIT. 


There is an error handling routine at program line 13000 in 
case someone presses control-C or gives a simple carriage return 
when asked for a nuaber. This simply prints an error message and 
returns to the main menu. The proper exit from the prograa is 
choice 8 on the wenu, "QUIT". That choice goes ta line 20000 af 
the progras, which does a "RESET" and an "FP" to escape the high 
resolution character generator and get the Apple back isto a 
standard configuration. The user should type in RUN MENU at this 
point to get a new progras. 
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IDEAL GAS CYCLES grew too big for memory. The first thing 
that was discarded was several REM’s. Fortunately, the internal 
program notation is similar to a physics text 50 a person familiar 
with thermodynamics can follow it. The second thing discarded was 
a lot of the text messages, such as instructions. These are out on 
disk storage until the prograa needs thea. This program uses all 
of the aystericus binary and text files on the thermodynamics 
disk, 


MOLECULAR MOTION AND GAS PRESSURE is another program that uses 
tricks that could potentially create trouble. The program does 
most of its display on the HGR2 screen. It draws text on this 
screen with shape tables. If there is a control-C€ or a request to 
"8" (quit) then the progras jumps to line 280 which resets TEXT and 
stops the prograa. 


The program has two sachine language portions appended to its 
end. The first is a shape table used to draw text on the high 
resolution screen, 


The second is a machine language subroutine called "MOVE" in 
line 100. This does several things very quickly. (1) It aoves all 
NA molecules ahead one time step, calculating the new positions and 
plotting thes. (2) It increments the instantaneous pressure PZ by 
adding to it the speed of every particle that hits the right wall 
of the cylinder. (3) It increments the tiae TI2 and shifts the 
graph left aif the time is greater then 248. This subroutine 
INSISTS that the variables in line 1000 be defined first before any 
other variables in aernory. The "MOVE" routine needs to find these 
variables in a certain place. 


Positions (X%, YZ) and speeds (V%, W%) have been scaled to 
maintain precision, They are plotted 14 times smaller than the way 
they are stored in memory. That is, X = ¥ + V for moving particles 
precisely but HPLOT 4/14, Y/14 is done to display the particles. 


MENORY USAGE: 2048 - 10000 (roughly) = main Basic prograe 
16384- 24576 = HBR2 screen for the display 
24576- 26199 = machine language routine to 
move all the sclecules 
26200- 31000 (about) = variable storage 
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SUGBESTED CHANGES 


The pressure is smoothed 
out in this line: 


HQ PS = (9 € Ph + PZ)? tO 


PZ is the pressure from molecules that hit the pressure gauge 
recently. PS is a longer time averaged pressure. If line 110 is 
deleted then the pressure graph will have acre spikes on it, On the 
other hand, the amount of time averaging can be increased: for a 
sacother graph. For example, PS = (P5 # 20+ PY) / 21 will make a 
very smooth pressure graph because it gives little weight to the 
few pressure fluctuation Px. (The molecule moving subroutine 
"MOVE" internally s@oothes out the pressure PY by taking the hits 
from molecules plus {5 times the old P% and dividing by 14.) 


The command keys are in lines 205-225 and 800-820, They can he 
other letters, 


The factor T2 in lines 620 and 450 controls how auch the 
teaperature changes every time a key is pressed. It can he bigger 
for faster changes. 


The gaxiaus number of particles is 127, The MOVE routine can’t 
handle any sore. 


The mimimum temperature is 2 degrees in line 450, If T is 
below 2 then the program won’t let the gas get colder, This is 
because it has trouble restarting the aolecules after they have 
stopped at absolute zero. The maximum temperature is 2500 degrees 
in line 400. There is no good reason for this particular number, 
except that at higher temperatures the molecules sometimes bounce 
from top to hottom of the cylinder in a single tiae step and 
confuse the MOVE routine. 
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PROGRAM STRUCTURE - MOLECULAR MOTION AND GAS PRESSURE 
Lines 10-20 draw the text string M€ on the high resolution 
screen at cursor position CH, CV 
Lines 100-190 are the main moving and plotting loop. 
‘The program spends aost of its tiae in here 


letting the gas sove and waiting for a user command, 


Lines 200-225 pick the place to go after a user presses 
one of the allowed commands. 


Lines 300-340 add one more aolecule at the center of the 
gas cylinder, 


Lines 344-348 clean up the screen after an error gessaqe, 
Lines 350-395 regove one particle, the last one nuebered NZ, 


Lines 400-440 calculate and print the temperature. They are 
only used once when the progras is initiallizing. 


Lines 600-640 aake the gas hotter. They multiply the speed of 
every particle by the factor 72. They multiply the 
20 segments of the Maxwellian distribution by 72, 
They also multiply the teaperature by T2 squared. 


Lines 650-470 sake the gas colder by the factor 12, 
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PROGRAM STRUCTURE - MOLECULAR MOTION AND GAS PRESSURE 


Lines 700-720 set up the nuaber bins in the 
Maxwellian distribution. 


Lines 750-770 set the speed and direction of a new particle (J) 
Lines 800-840 draw the command list on the screen. 


Lines 1000-3200 print the title, draw the display, 
and initiallize al] of the constants. 


Line 43999 sets up a pointer to the ASCII shape table that is 
appended to the end of line 63999, 


63999 79 = PEEK (121) + 256 # PEEK (122):29 = 79 + 73: 
POKE 232, 29 - INT (29/254) & 256: 
POKE 233, 79 / 256 : RETURN :REM SET UP SHAPE TARLE 


- - ~ HIDDEN SHAPE TABLE FOLLOWS AFTER THIS LINE - - - 
MOSTLY ASCIT LETTERS FOR HI-RESOLUTION TEXTS 


- - ~ A MACHINE LANGUAGE PARTICLE MOVING ROUTINE - - - 
FOLLOWS AFTER THE SHAPE TABLE 


A scene from Molecular Motion and Gas Pressure: 
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COMMANDS : iy 
L = Less M ES 
i = HOTTER "€"= COLDER 
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